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The possibility of explaining the positron and electron excess recently found by the PAMELA 
and ATIC collaborations in terms of dark matter (DM) annihilation has attracted considerable 
attention. Models surviving bounds from, e.g, antiproton production generally fall into two classes, 
where either DM annihilates directly with a large branching fraction into light leptons, or, as in the 
recent models of Arkani-Hamed et al., and of Nomura and Thaler, the annihilation gives low-mass 
(pseudo)scalars or vectors cj} which then decay into fi'^ n~ or e'^e~ . While the constraints on the first 
kind of models have recently been treated by several authors, we study here specifically models of 
the second type which rely on an efficient Sommerfeld enhancement in order to obtain the necessary 
boost in the annihilation cross section. We compute the photon flux generated by QED radiative 
corrections to the decay of and show that this indeed gives a rather spectacular broad peak in 
E^da/dE, which for these extreme values of the cross section violates gamma-ray observations of 
the Galactic center for DM density profiles steeper than that of Navarro, Frenk and White. The 
most stringent constraint comes from the comparison of the predicted synchrotron radiation in the 
central part of the Galaxy with radio observations of Sgr A*. For the most commonly adopted 
DM profiles, the models that provide a good fit to the PAMELA and ATIC data are ruled out, 
unless there are physical processes that boost the local anti-matter fluxes more than one order of 
magnitude, while not affecting the gamma-ray or radio fluxes. 



There have recently been indications of a very inter- 
esting enhancement in the amount of cosmic ray elec- 
trons and positrons detected near the Earth, both seen 
by PAMELA in the ratio of positrons to the sum of elec- 
trons and positrons between a few GeV and 100 GeV 
Q, and by ATIC in the sum of electrons and positrons 
at several hundred GeV to 1 TeV While these so 
far unexplained excesses might be due to standard as- 
trophysical processes Q , positrons also constitute one of 
the promising channels in which to search for dark matter 
(DM; for reviews, see Q), and these new experimental 
findings have therefore already triggered a large number 
of theoretical analyses trying to explain the data as being 
induced by DM annihilation or decay (see e.g. Ref. 5| 
and references therein for supersymmetric DM, Refs. [6| 
for alternative DM scenarios and Refs. Q for decaying 
DM scenarios). In general, these analyses seem to point 
at the need for DM particles with masses in the TeV 
range that annihilate, with a very large rate, dominantly 
into charged light leptons. 

The bremsstrahlung process, falling like E~^, is gener- 
ally regarded in particle physics as having a "soft" spec- 
trum. In the astrophysical context, this is, however, 
on the contrary a quite hard spectrum, since most of 
the background 7-ray spectra like those from accelera- 
tion near supernova remnants usually fall like E~'^ or 
faster. Gamma-rays from DM generally feature a spec- 
trum that is somewhere in between these two at low ener- 
gies {E~^-^) and drops even faster close to the DM parti- 
cle mass (for important exceptions see, however, [§].) 



If the DM particles x annihilate directly into a pair of 
charged leptons, the photon distribution from the process 
XX ~^ £^i~^, for rrix ^ rni, is to a good a ppr oximation 
of the Weizsacker- Williams form (see, e.g., jlOj): 
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where x = /m-^ and {av)u is the annihilation rate for 
the lowest order process XX ~^ i~^£~ (Note that the above 
approximation also breaks down when there is a symme- 
try that suppresses the annihilation into two-body, but 
not into three-body final states [U). 

This case has recently been treated by [13, [3, fl^ - 
(The last of these references also briefiy treats, but leaves 
for a more detailed calculation, the kind of processes we 
will compute here.) It was found that the gamma-rays 
produced in DM models with these annihilation modes 
lead to rather severe constraints. Even more stringent 
bounds on this type of DM models that try to explain 
the PAMELA and ATIC data arise from the synchrotron 
radiation produced by the resulting population of elec- 
trons and positrons, in realistic models of the DM density 
distribution and for a wide variety of assurnptions about 
the magnetic field in the inner Galaxy [l^, [l^ . 

It remains to consider another possibility, where DM 
annihilates into a new type of light (sub-GeV) particles 
^that in turn dominantly decay into light leptons (see 
[161 1 for a general account of this idea). The advantage 
of this type of models is that the strongly constrained 



2 



Arkani-Hamed et al. type 


Nomura- Thaler type 
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ms[GeV] 


ma[GeV] 


AHl 


0.1 


scalar 


100% 




iVl 


5 


0.5 


AH2 


0.1 


vector 


100% 




N2 


20 


0.36 


AH3 


0.25 


vector 


67% 


33% 


m 


20 


0.5 


AH4 


0.25 


scalar 




100% 


iV4 


20 


0.8 












iV5 


50 


0.5 



TABLE I: Our benchmark scenarios. 



decay into hadronic modes (see, e.g., [17|) is kincmat- 
ically forbidden and that Sommerfeld enhancements in 
the hmit of the small galactic DM velocities expected to- 
day allow for the very large annihilation cross sections 
that are needed to explain the PAMELA/ATIC results, 
but which at first seem to be at odds with the cross 
sections required to get the right thermal relic density 
for the DM. Another interesting feature of the Arkani- 
Hamed et al. model [l^ is that it encompasses ideas that 
have been proposed to explain the WMAP haze 
the INTEGRAL excess [H]. As pointed out in 



18 1 and 
16,13, 



one may basically distinguish between scalar and vector 
(j) and whether or not < 2m^ (in which case it dom- 
inantly decays into e^e^). For > m7r, even decays 
into pions should be taken into account (which we ne- 
glect here). While > lOMeV is roughly needed not 
to be in conflict with Big Bang Nuclesynthesis, one has 
to require > 100 MeV in order to get Sommerfeld 
enhancements of the order 10"^ — 10"* that are needed to 
explain the PAMELA/ATIC result with these types of 
DM models. Based on this discussion, we adopt the four 
benchmark settings A1-A4 summarized in Tab. [J 

While [3l describes a rather general set-up, [2l[ in- 
troduces a concrete realization of this idea; the proposed 
model has the appealing feature of containing a "stan- 
dard" Peccei-Quinn axion and can be embedded in a 
fully realistic supersymmetric scenario. Here, DM an- 
nihilates into a scalar s and a pseudoscalar a, XX ~^ sa. 
With a mass scale of 360 MeV < < 800 MeV, the lat- 
ter mostly decays into muons, which subsequently decay 
into electrons or positrons. The benchmark models for 
this setup N1-N5 are also given in Tab. HI 

For the first a particle created in the xX annihila- 
tion, we analytically compute the photon multiplicity 
{dN / dE^)^"'^ from a ^ /i"*"^ 7 in the rest frame of a. 
We then make a Lorentz boost back to the DM frame, 
i.e. the Galactic rest frame, to get 
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with 7 = {m^/ma) [l — {ml — m^)/(4m^)] since the an- 
nihilation takes place essentially at rest (typical galactic 
velocities are 10~^). Axions resulting from s — s- aa we 
treat in a similar way, boosting them first to the s-frame 
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FIG. 1; The various possible photon spectra that can arise 
from DM annihilating to new light particles which in turn 
decay into charged leptons. For the models A''! - A'^5, we 
neglect here the decay of s to tau-leptons or bottom quarks - 
see Fig.[5]for an example of how this changes the spectra. For 
comparison, we also indicate the spectrum from DM directly 
annihilating to charged leptons. 



and from this to the DM frame. Since s may have a 
mass up to 50 GeV, the gamma-ray spectrum may even 
receive important contributions from its decay into bot- 
tom quarks or tau leptons, a possibility which we will 
shortly return to. (Bremsstrahlung from electrons in the 
muon decay will give 7s of lower energies and will thus 
not be important for our constraints.) 

Summing up all these contributions, we arrive at the 
total photon spectrum in the DM frame that we show 
in Fig. [2 for the models N1-N5 in Tab. HI We also in- 
clude the corresponding spectra obtained in the Arkani- 
Hamed et al. set-up (models A1-A4) and, for comparison, 
the case of 1 TeV DM particles directly annihilating into 
e+e~ or /i^/i^. Please note that, from Eq. ([2|), the quan- 
tity dN/dx for the models listed in Tab. U is independent 
of as long as 3> rua, ms\ the direct annihilation of 
DM into leptons, on the other hand, does contain a log- 
arithmic dependence on m^. Let us mention that while 
Eq. ([1]) provides a rather good approximation to our an- 
alytic results for photons radiated from e^e~ pairs, it 
overestimates the photon yield from muons (especially 
when the mass of the decaying particle is close to 
like, e.g., in model AH4). 

Once a DM profile p{r) is assumed, it is straightforward 
to estimate the corresponding gamma-ray flux from a 
solid angle Afi towards the galactic center: 
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where A is the line of sight distance. In Fig. [21 we 
compare the resulting flux to the gamma-ray data from 
the galactic center taken by the H.E.S.S. telescope [2^ . 




FIG. 2: The total gamma-ray spectrum dN/dE^, for an NFW 
halo, from a 1 TeV DM particle annihilating into a pseu- 
doscalar a (decaying to muons) and a scalar s which decays 
to aa (solid line) or only in 95% of the cases into aa and in 
5% into bb (dotted line) or t'^t^ (dashed line). The masses 
for a and s are those of model N3 of Tab. |IJ so the solid line 
corresponds to the A'^S-line shown in Fig. [T] 



FIG. 3: Exclusion plot in the av vs. mass plane. The two 
sets of curves give the maximum annihilation cross section 
compatible with radio observations of Sgr A* for Einasto and 
NFW profiles. The color code of the curves is the same as in 
Fig.m The shaded region, corresponding to the range of anni- 
hilation cross sections that provide a good fit to the PAMELA 
and ATIC data, appears to be in conflict with observations, 
unless the DM profile is more shallow than Einasto. 



which has an angular resolution of about 0.1°, thus 
Aft = 10~^ sr. We here show the spectrum for model 
N3 and, for comparison, the case wliere s decays not 
only to axions but with a branching ratio of 5% to bb or 
T~^T~ (which is the typical case for the model presented in 
PH). By comparison with Fig. [1] it is straightforward to 
arrive at the corresponding spectra for the other models 
in Tab. HI We have here adopted a so-called Navarro- 
Frenk- White profile [2^ , with the same parameters as in 
Ref. [l3|. Note that the gamma-ray spectra in this case 
are consistent with the HESS data, unlike the case of the 



annihilation modes discussed in [l3j , for the same density 
profile. Assuming a profile p{r) oc r~^''^, as needed to ex- 
plain the WMAP 'Haze' (see Ref. [ll|), the constraints 
become much more stringent. However, at the same time 
they become much more sensitive to the dependence of 
av on the velocity dispersion of DM, which inevitably 
increases in the vicinity of the supermassive black hole 
at the Galactic center. As we shall see soon, however, 
it is possible to derive even tighter constraints without 
making assumptions on the small- v behaviour of av. 

Before that, however, let us note that another potential 
source of gamma rays from DM annihilations are dwarf 
galaxies, like the Sagittarius dwarf galaxy, observed by 
HESS [il]. The HESS observations put an upper bound 
on the integrated gamma flux above 250 GeV of < 
3.6 X 10^^^ cm^^s^^. Assuming an NFW (isothermal) 
profile in the Sagittarius dwarf galaxy, this can be trans- 
lated to the limit av < 7.4 x lO"^^ (2.2 x IQ-^S) cm^ s^^ 
for model N3. For the other models in Tab. [H the lim- 
its differ by a factor of a few as indicted by the spectra 
in Fig. [T] For other dwarf galaxies, the limits are sim- 



ilar: using a conservative estimate of the line of sight 
integral from Ref. psj . the limits on the gamma flux 
from Willman 1 as observed by Magic [1^, e.g., translate 
to av < 1.3 x 10^^^ cm'^ s~^. However, the uncertain- 
tics from dynamical constraints |25j are large and im- 
proved future data might result in better constraints. As 
one typically needs a boost of order 10'^ to explain the 
PAMELA data, we note that the limits derived here are 
very close to the required av. This means that for some 
models, like AH1-AH3, the more optimistic scenarios for 
the halo profile of e.g. the Sagittarius dwarf are excluded. 

A rather stringent constraint on the rate of injection 
of high energy in the Galaxy comes from the analysis 
of the synchrotron radiation produced by these particles 
as they propagate in the Galactic magnetic field. Al- 
though observations of different targets and at different 
wavelengths provide interesting constraints [27| , the most 
stringent ones come from radio observations of the Galac- 
tic center, where the DM density is highest 1^, 27, 2^. 

The synchrotron luminosity generated by a distribu- 
tion of electrons and positrons produced by a DM distri- 
bution with profile p{r) in a magnetic field B{r) is 



vLiy = 2TT^r I drr^ p\r)E„YJE„) 



(4) 



where Ep = iy^^^[0.29{3/4:7r){e/meC^fB{r)]-^/^, 
Ye{E) = f^"" dE' dNe/dE' and we have adopted the 
monochromatic approximation for the synchrotron 
emission, assuming P{i',E) = (87r/9-\/3) J(i^/z^c — 0.29), 
with Vc = {3eB E'^) / (ATTmlc^) , for its spectrum. 

By comparing the predicted synchrotron radiation 
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with radio observations, we can set limits on the anni- 
hilation cross section for any given annihilation channel, 
following a procedure similar to Ref. ■ The most strin- 
gent constraint comes from the upper limit on the radio 
emission from a cone with half-aperture of 4" towards Sgr 
A* at = 0.408 GHz which we translate in Fig.[2]to 
the av vs. mass plane. Let us stress that the av plotted in 
Fig. [3] is the effective annihilation cross section, including 
both Sommerfeld enhancements and boosts due to sub- 
structures. The only way to avoid our constraints would 
thus be to boost the local anti-matter fluxes by more than 
one order of magnitude without affecting the gamma-ray 
or radio fluxes. Although this theoretical possibility can- 
not be ruled out (e.g. Refs.[30l|). it appears to be unlikely 
for a realistic distribution of substructures in the Milky 
Way halo. Numerical simulations seem to indicate that 
the boost factors due to substructure is rather small [3l[ . 
How big the boost factors could be are still under debate 
and recent simulations [3l| indicate that locally they are 
at most a factor of a few. A recent study [sll develops a 
model that indicates that the local boost could be about 
a factor of ten. The details of the mechanism giving such 
large boosts arc yet to be presented, however. For more 
discussion about boost factors, see Ref. [l3l |. 

The two sets of curves give the maximum annihila- 
tion cross section compatible with radio observations of 
Sgr A* for two different DM profiles: Einasto and NFW. 
The shaded region, corresponding to the range of an- 
nihilation cross sections that provide a good fit to the 
PAMELA and ATIC data, appears to be in conflict with 
observations, unless the DM profile is more shallow than 
expected in current models of structure formation. How- 
ever, if the DM interpretation of the PAMELA data was 
corroborated by additional evidence, then our result can 
be interpreted as a hint of the shallowness of the DM 
profile. 

Profiles steeper than NFW - like the p(r) cx r"^-^ 
needed to explain the WMAP 'Haze' ^ are ruled 
out by a rather larger margin. This confirms the dra- 
matic importance of the multi-wavelength approach to 
DM studies [H, [H, [13, El], especially for DM models 



tailored to explain anomalies in astrophysical observa- 
tions. 
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